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Целью работы является создание модели сложного распреде- 
ленного объекта — цементного производства на основе 
раскрашенных временных сетей Петри, ее анализ и оценка с 
точки зрения эффективности и достоверности. В работе 
рассмотрены различные виды сетей Петри, способы их 
задания и динамика работы. Для раскрашенных временных 
сетей Петри сформулирован алгоритм функционирования в 
матричной форме. На основе раскрашенных временных сетей 
построена математическая модель, позволяющая спрогнози- 
ровать объем цементного производства за определенный 
период. Модель реализована на языке программирования 
С++. Проведено сравнение результатов моделирования с 
фактическими данными. Установлено, что реализованная 
модель с достаточной точностью предсказывает объем 
выпуска продукции цементного производства. Разработанные 
методы работоспособны и применимы в моделировании 
производственных процессов в составе автоматизированной 
системы управления технологическими показателями. 
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Введение. В работе рассмотрена методика создания модели на примере сложной распределенной системы — 


технологического процесса функционирования цементного производства. Процесс выпуска цемента является одним из 


примеров недетерминированных динамических параллельных производственных систем, проблема моделирования 


которых связана как с возможной хаотичностью системы, так ис необходимостью учитывать динамику подсистем. 


Для описания и анализа таких систем могут применяться сети Петри [1-3] и их разновидности, например, не- 


четкие [4], временные [5—6], раскрашенные [7]. 


В [8-10] рассмотрены окрестностные модели обжига клинкера цементного производства. В [11] при модели- 


ровании процесса функционирования цементного производства использованы временные сети Петри, достоинствами 


которых являются динамическое отражение состояний моделируемой системы и возможность анализа свойств полу- 


ченной модели. 


* Работа выполнена при финансовой поддержке РФФИ (проект №16-07-00-854). 
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В работе рассмотрены раскрашенные временные сети Петри, являющиеся сетями более высокого уровня и 
позволяющие, по сравнению с обычными сетями Петри, анализировать дополнительные свойства моделируемых 
процессов без усложнения структуры сети. 

Так, во временных сетях, в отличие от обычных сетей Петри, переходы срабатывают с некоторой задержкой, а 
маркеры находятся в позициях определенное время, что дает возможность моделирования не только последовательно- 
сти событий, но и их привязку ко времени. 

Раскрашенные сети Петри позволяют одновременно моделировать несколько параллельных потоков различ- 
ных материалов или событий в процессе функционирования сложных систем. В аналогичных моделях на основе 
обычных сетей Петри приходится искусственно вводить дополнительные позиции, не являющиеся отображениями 
элементов процесса, служащие для упорядочения запусков переходов сети и разделения материалов или событий, что 
усложняет пространственную структуру модели и затрудняет ее интерпретацию. 

Раскрашенные временные сети Петри, используемые в работе, объединяют в себе приведенные достоинства 
как раскрашенных, так и временных сетей. 

В работе построен опытный образец модели расчёта объема выпуска продукции цементного производства по 
месяцам в течение одного года на основе раскрашенных временных сетей Петри. 

Разработана программа на языке С--, позволяющая рассчитать производительность цементного производства 
за заданный период. Реализованная модель с достаточной точностью предсказывает объем выпуска продукции. 

Способы задания и правила функционирования сетей Петри. Существует три эквивалентных способа за- 
дания сети Петри: графический, аналитический и матричный [3, 11]. 

Графически сети Петри представляются в виде двудольных графов. Множество вершин состоит из непересе- 


кающихся подмножеств позиций Р={р,}, 1=1....,п и переходов Т = {1,}, /=1....,т, а множество дуг разделяется на 
два подмножества {(р;,1 5} сРхХТ и{(1 _ р)}<ТхР. В изображении графов, представляющих сети Петри, позиции 


обозначаются кружками, а переходы — планками. 
Далее рассмотрим аналитическо-матричный способ задания сетей Петри [11]. Сеть Петри задается следую- 


щим набором РМ =(Р,Т,К`_,В*, Но), где: 

-Р={р,ро,....Р„} — конечное непустое множество позиций; 

- Т={1,6,...м} — конечное непустое множество переходов (множества Р и Т не пересекаются: 
РПТ=@); 

— К еВ"”" — матрица инцидентности дуг, входящих в переходы; 


- ®*еВ"”" — матрица инцидентности дуг, выходящих из переходов; 


— Ш= {ш, Пао} — вектор начальной маркировки сети Петри. 


Приведем алгоритм функционирования сети Петри: 
1. Текущая маркировка сети равна начальной и = о. 


2. Переход 1, (]=1...т) при текущей маркировке и разрешен, если п2е,.К_, где 


г; =[0,0,...,1:;..Ое В” — строка, содержащая нули везде, за исключением ]-го элемента. Заметим, что разрешен- 


ных переходов при текущей маркировке может быть несколько. Если нет разрешенных переходов, сеть достигла 
тупиковой маркировки, дальнейшее функционирование невозможно, конец алгоритма. Иначе переходим к пункту 3. 


3. Случайным образом выбирается один из разрешенных переходов {, (1=1....,т). 
4. Маркеры перемещаются из входных позиций выбранного перехода 1, во все его выходные позиции по 
формуле и=и+ е; .Ю,где =А*-В — матрица инцидентности сети Петри. 


Виды сетей Петри. Раскрашенная временная сеть Петри. Во временных сетях Петри [5] 


РМ = (РТ, К, В", щ,7,5) вводятся в рассмотрение временные задержки маркеров в позициях и время срабатыва- 


ния разрешенных переходов, где: 


-5={51,5,,...,5„} — вектор задержек маркеров в позициях; 


-П= бойне а р — вектор времени срабатывания разрешенных переходов. 
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Раскрашенная сеть Петри [7] РМ. = (Р,Т, С.К, К, йо) отличается от РМ наличием цветов, матрицей, а не 
вектором, начальной маркировки и блочной структурой матриц инцидентности, где: 


— С={<1,со,.... са} — цвета сети; 


_ В еВ“"" блочная матрица инцидентности дуг, входящих в переходы; 


В = ‚ где К, — матрица инцидентности входящих дуг цвета С’, К=1...., а; 


_-В*е В“ блочная матрица инцидентности дуг, выходящих из переходов; 


В" =| 2 | где в, — матрица инцидентности выходящих дуг цвета С, ; 








. 
[Ки | 
—Й0еВ“”” — матрица начальной маркировки. 


Обобщением временной и раскрашенной сети является раскрашенная временная сеть Петри 
р а. 
РМ. =(Р.Т,С,К ‚К, Во,2,5). 


Представление цементного производства посредством раскрашенных временных сетей Петри. Рас- 
смотрим в данном пункте реализацию модели на примере сложного распределенного объекта — технологического 
процесса функционирования цементного производства ЗАО «Липецкцемент». 

На рис. | изображен граф раскрашенной временной сети Петри, иллюстрирующий производственный цикл 
цементного производства. 





Рис. 1. Граф раскрашенной временной сети Петри цементного производства 


Позиции Ру- Р5 на рис. 1 соответствуют следующим складам: 


— ри — склад сырья; 


р› — силос сырьевой муки; 

— рз — силос сырьевой муки; 

— р. — склад клинкера; 

— р5 — силос цемента. 

Переходы 1-15 соответствуют агрегатам: 


— 1-14 — сепараторные мельницы 3,2 Х 8,5 м; 
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1; —& — трубные мельницы 4,2 Х 10 м; 

— 1-4 — вращающиеся печи 4Х 60 м; 

— 1, — вращающаяся печь 5 Х 75 м; 

— по-йз3 — цементные мельницы 3 Х 14 м; 

— 14-6 — Цементные мельницы 3,2Х 15 м. 

Граф сети Петри имеет цвета С = {С;, Съ, Су, Сл, С5, С} (рис. 1). 


Блочные матрицы инцидентности входов и выходов переходов К и К* соответственно равны: 








К К 
К = ы. .®= ® 
К К 
В частности, для цвета С› матрица А› равна: 
о о0бо0о00 0 0 000000001 
оо0о000=я хх 00000000 
=ооо0о000 0 0 00000000, 
Фоо0о000 о 0 00000000 
[001 соб бох об обо сбооо 
где г. —Щ количество сырья, поступающего со склада Р› во вращающуюся печь ГА 7 ( 7 ='7, 8) . Остальные элементы 


матрицы А равны нулю, так как в графе сети Петри имеется только две входящие оранжевые дуги С› ‚ соответству- 
ющие переходам 1, —& и позиции ро (рис. 1). 
Матрица > равна: 


ооо о 


2+ 2+ 2+ 2+ 
р 72 РЗ 14 


сое >= 
всеее >= 
веоее = 
сое >= 
всеее >= 
сое >= 
всеее = 
всеоее >= 
веоеозе >= 
всее >= 
сеое >= 
веоее >= 








№Ю=о о о 0 
0 0 0 0 
[о 0 0 0 .| 
где ИН — производительность сепараторной мельницы 1, (]=1....,4). Другие элементы матрицы Е> являются 


нулевыми, так как в графе имеется только четыре выходящие оранжевые дуги С. , соответствующие переходам 1- 14 
и позиции р» (рис. 1). 
Остальные матрицы Кг, К; (К=1,...,6) формируются аналогично. 


Раскрашенная временная сеть Петри цементного производства функционирует по следующему алгоритму: 
1. Начальное время функционирования сети т=0 ‚; начальная маркировка цу описывает количество материала 


на складах в начальный момент времени; время функционирования производства равно Т часов. Все переходы сети 
Петри ЯВЛЯЮТСЯ незаблокированными. 


2. Незаблокированные переходы сети {,, ]=1,..,т’ последовательно проверяются на разрешенность. Пере- 
ход 1, при текущей маркировке и, разрешен, если и. 2е,®К’,гдее, =[0,0,...,1,,...,0] — строка, содержащая нули 


везде, за исключением го элемента; операция 69 означает произведение строки е; на каждую из матриц К ь обра- 


Т 
зующих блочную матрицу К = [к ‚К Я . Если нет разрешенных переходов, переход к п.4. 
3. Маркеры перемещаются в разрешенный переход 1 г Результат начала запуска перехода 1 ; при текущей 


маркировке и. записывается как и. =ИК. —е; 9. Далее переход блокируется на время э; выполнения операции. 
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Переход к п.2. 
4. Сдвиг времени на т=т+1 час. Если т>Т, то алгоритм завершен. Маркеры переходят в выходные позиции 


разрешенных в п.3 незаблокированных в данный момент времени переходов {, по следующей формуле 


Шен. +е, ® К’. Переход к п.2. 


По данному алгоритму была разработана программа на языке С++, целью которой является моделирование 
динамики состояний системы и вычисление объема выпуска продукции цементного производства за определенный 
период. 

Рассмотрено функционирование модели на основе реальных данных за 2012 год. Результаты, полученные в 
процессе моделирования, приведены в табл. 1. 





Таблица 1 
Относительная ошибка выпуска продукции 
а Относительная ошибка 
моделирования, % 

январь 0,134940 
февраль 0,480579 
март 0,691856 
апрель 0.080636 
май 0,257613 
ИЮНЬ 0.099582 
ИЮЛЬ 0,235929 
август 0,036875 
сентябрь 0,076327 
октябрь 0,407839 
ноябрь 0,452516 
декабрь 0,071766 











Относительная ошибка найденного объема выпуска продукции вычисляется по формуле: 


бР, = УТУ *100% , 
У 
где у; — реальные данные выпуска продукции за 1-ый месяц (= 1,...,12); 
у; — модельные данные выпуска продукции за 1 -ый месяц. 


Средняя относительная ошибка моделирования составляет 0,25% и является допустимой для применения 
предложенной модели при прогнозировании объема выпуска продукции цементного производства. Таким образом, 
проведенные расчеты свидетельствуют об адекватности разработанной модели. 


Заключение. В работе реализовано представление раскрашенных временных сетей Петри в матричной форме, 
сформулирован алгоритм их функционирования. 

Раскрашенные временные сети Петри применены в моделировании сложного распределенного объекта — це- 
ментного производства. В качестве приложения рассматривается ЗАО «Липецкцемент». 

Разработано программное обеспечение на языке С++ для реализации модели функционирования цементного 
производства. Проведено сравнение результатов моделирования с фактическими данными. Оценена пригодность 
разработанной модели для прогнозирования производственных процессов на основе данных 2012 года. 

Реализованная модель с достаточной точностью предсказывает объем выпуска продукции и может быть эф- 
фективно использована для прогноза и анализа динамики производственных процессов. 


Библиографический список 

1. Питерсон, Дж. Теория сетей Петри и моделирование систем / Дж Питерсон. — Москва : Мир, 1984. — 264 с. 

2. Котов, В. Е. Сети Петри / В. Е. Котов. — Москва : Наука, 1984. — 160 с. 

3. Васильев, В. В. Сети Петри: параллельные алгоритмы и модели мультипроцессорных систем / В.В. Васи- 
льев, В. В. Кузьмук. — Киев : Наукова думка, 1990. — 213 с. 


Седых И. А. и др. Применение раскрашенных временных сетей Петри 





4. Леоненков, А. В. Нечеткое моделирование в среде МАТГАВ и Ё22уТЕСН / А. В. Леоненков. — Санкт Пе- 
тербург : БХВ-Петербург, 2005. — 736 с. 

5. Воевода, А. А. Временные сети Петри и диаграммы ОМГ / А. А. Воевода // Науч. вестн. НГТУ. — 2009. — 
№ 4(37). — С. 169-174. 

6. \\апь, 7. Типед Рей №5: ТБеогу ап АррИсаноп / Г. У’апз. — Могмей : Ка\ег Аса4епис РиБИзВегз, 1998. 
— 296р. 

7. Ехлаков Ю. П. Цветные сети Петри в моделировании социально-экономических систем / 
Ю. П. Ехлаков, В. Ф. Тарасенко, О. И. Жуковский // Доклады Том. гос. ун-та систем управления и радиоэлектроники. 
— 2013. — №3 (29). — С. 83-92. 

8. Шмырин, А. М. Наличие экстремумов параметрического уравнения печи обжига клинкера / А. М. Шмы- 
рин, И. А. Седых, А. П. Щербаков, А. Г. Ярцев // Вести высших учебных заведений Черноземья. — 2015. — №1(39). — 
С. 62-67. 

9. Шмырин, А. М. Исследование окрестностной модели печи обжига клинкера с учетом допустимых значе- 
ний параметров / А. М. Шмырин, И. А. Седых, А. П Щербаков, А. Г. Ярцев // Вестник Липецкого Государственного 
Технического университета. — 2015. — №2(24). — С. 11-14. 

10. Шмырин, А. М. Алгоритмы идентификации и управления функционированием окрестностных систем, по- 
лученных на основе сетей Петри / А. М. Шмырин, И. А. Седых // Управление большими системами. — 2009. — 
Вып. 24. — С. 18-33. 

11. Блюмин, С. Л. Окрестностное моделирование сетей Петри / С. Л. Блюмин, А. М. Шмырин, И. А. Седых, 
В. Ю. Филоненко. — Липецк : ЛЭГИ, 2010. — 124 с. 


ВегГегепсе5 

1. Раегзоп, 7. Теопуа зееу Рей1 1 шодейгоуаше зузет. [Рец Ме! Твеогу ап Фе МодеПпе оЁ Зу$етл$.] Мозсом: 
МИ, 1984, 264 р. ап Визяап). 

2. Котоу, У. Е. 5еи Ревч. [Рей пез.] Мозсо\: МаикКа, 1984, 160 р. (т Виззап). 

3. УазШеу, У. У., Килшак У. У. Зей Ребт: рагаПеГпуе а1опёту 1 тодей шаГаргозеззогиуКВ зузет. [Рейт пе: 
рагаПе] а12огИил$ ап4 по4е]$ оЁ ми@ргосеззог зубепаз.] Клеу: МачКоуа ЧитКа, 1990, 213 р. (т Визчап). 

4. Геопепкоу, А.У. Месве кое то4ейгоуаше у этеде МАТГАВ 1 в 27уУТЕСН. [Еа77у тодейпо ш МАТГАВ апа 
Ё177УТЕСН зой\аге епугоптепе.] 51.РеегзЬаго: ВКВУ-Раегого, 2005, 736 р. (п Кизз1ап). 

5. Уоуеуода, А. А. Угетеппуе зен Рейл 1 Ф1азтатту ОМГ. [Титед Реёл пе апа ОМГ, аснуцу Фаотатз.] Зс1еп- 
ибс ВиПейп оЁ МТО, 2009, по. 4(37), рр. 169-174 (ш Виззлап). 

6. \\!апо, 7. Титеа Рейл Ме{з: ТВеогу ап АррПсаНоп. Мог\меП: КТа\мег Асадепис РиБИзВегз, 1998, 296 р. 

7. ЕКШаКоу, У. Р., Тагазепко, У. Е., УБиКоузКу, О. 1. Тзуешуе зей Рейт у тодейгоуапи зоблаГпо- 
екопописвезЮКЬ з15{ет. [Со]ог Рей пеё5 ш шо4деНпе оЁ 5ос1ю-есопопис зууетз.] Ргосее 4115 о ТОЗОК Ушуегзиу, 2013, 
по. 3 (29), рр. 83-92 (т Виззап). 

8. оптуги, А. М., Зедуки, Г. А., ЭсВетфакКоу, А. Р., Уанбеу, А. О. МайсШе екзметитоу рагателевезКо2о 
игаупетуа рес об7Шра КНикКега. [Те ргезепсе оЁ сПпКег Рагпасе рагатейлс едиайоп ехшета.] Ме\из оЁ Н1эВег ЕдисаНопа| 
зи @оп$ ое Срегпо7ет Вестоп, 2015, по. 1(39), рр. 62—67 (т ВКиззап). 

9. эытуг, А. М., ЗеаукЬ, Т. А., ЗсНетРаКоу, А. Р., Уаивех, А. С. [3едоуаше оКгезпозтоу шодей ресШ об7Шга 
КПоКега $ асбеют доризитуКЬ 7пасбешу рагатешоу. [КезеагсВ оЁ а пеебогроо4 плоде! оЁ а сПпКег КПп факто шю 
ассоип( ап; 1Ые рагатеег уашез.] Уезиик ГЗТО, 2015, по. 2(24), рр. 11-14 (т Виззап). 

10. Зытуни, А. М., ЗедукВ Т. А. АШгомыту 14епННКаи 1 ирга/ешуа Рай 1юоптоуатет оКтезепозетукН з15ет, 
росвеппукИ па озпоуе зеёеу Рейт. [А]е отт оЁ1Аепийсайоп апа орегайопа! сопйго] оЁ пеоНботтоо4 зузетз БийЕ оп Фе 
Баз1$ оЁ Рец пе.] Гагое-зсае Зу\епа$ Сопго|, 2009, 155. 24, рр. 18 — 33 (ш Ваззап). 

11. ВШуипит, 5.Г., ЗЬтупа, А.М., ЗедукКВ, Т.А., ЕПопепКо, У.У. ОКгезтозтое то4еПгоуаше зе&еу Рец. [ПМе1эВ- 
Богпоо4 то4деНпз оЁ Рен пез.] Глрезк: ГЕСТ, 2010, 124 р. (т Ваззлап). 


Поступила в редакцию 14.07.2016 Кесетуеа 14.07.2016 
Сдана в редакцию 14.07.2016 бибтишеа 14.07.2016 
Запланирована в номер 30.09.2016 бсведще4 1 Ше 1550е 30.09.2016 


